Appln. No. 09/745.590 

Response dated j£ » f . 2004 

Reply to Office Action of September 8. 2004 


REMARKS/ARGUMENTS 

Claims 1-36 are pending. Claims 1-31 have been rejected, but 
Applicants have presumed that Claims 1-36 were meant to be rejected in the Final 
Rejection dated September 8 ? 2004. New claim 37. which is dependent from Claim 1. 
has been added and is supported on page 8, lines 3-7. Claims 32 and 33 have been 
amended to cancel an inadvertent "t" present in these claims. 

THE §112 REJECTION 

Claims 1-31 have been rejected under 35 U.S.C §112, second paragraph, as 
being indefinite for failing to point out and distinctly claim the subject matter which 
applicants regard as the invention. Applicants respectfully disagree and respond as 
follows. 

The Office Action gives two reasons as to why the claims are indefinite. The 
first reason as stated by the Office Action is: 

"Where applicant acts has his own lexicographer to 
specifically define a term of a claim contrary to its ordinary 
meaning, the written description must clearly redefine the 
claim term ... It appears that the 'box counting dimension' 
in the claims is used by the claim to mean 'packing 
density'. . . . The term is indefinite because the 
specification does not clearly redefine the term." 

The second reason as stated by the Office is: 
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"It appears from the reading of the specification that this 
[box counting dimension] is based on a software product 
and entails parameters of pore size and number of spheres, 
that is a pore channel size. The use of a non-standard, such 
as a software program is viewed as indefinite, because these 
programs may change over time, ..." 

Applicants use the box counting dimension of the surface structure as it is 
known in the art. That is, the box counting dimension has not been used in any way 
contrary to its ordinary meaning. The "box counting dimension" or also commonly 
referred to as the "fractal dimension" as determined by the box counting method of a 
surface is a well known parameter in heterogeneous catalysis. This is shown by the 
more than 200 citations given by a search of the open literature, which is attached. 
Typically, in heterogeneous catalysis the fractal dimension is determined on a scale 
much less than 1 micrometer by employing adsorbed gas molecules of different sizes, 
which is illustrated in the enclosed articles. Colloids and Surfaces, Physiochemical 
and Engineering Aspects, 187-188 (2001 ) 257-265 and Catalysis Letters, 2(1989) 
191-200. In addition, the box counting dimension or fractal dimension has been used 
and described in the patent art such as in U.S. Pat. Nos. 5.307,292; 5.444,795; 
5,666,169; 5,719,784; 5,848,177; 5,982,386 and 6,143,158. The box counting method 
is specifically described as using one known in the fractal art "such as those described 
by Falconer, Fractal Geometry Mathematical Foundations and Applications, John 
Wiley and Sons, NY," in which the pertinent section is enclosed, (page 6, lines 2-6 of 
the specification). 
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The only difference with the above cited heterogeneous catalysis art utilizing a 
fractal/box counting dimension is that, applicants have discovered a surprising result 
for the catalytic treatment of exhaust gas. The surprising result arises when the 
fractal/box counting dimension is at least 1 .2 for a catalyst-substrate surface wherein 
the fractal/box counting dimension is determined on a scale greater than 1 micrometer 
to about 1 mm. The box counting methodology does not change from the above cited 
heterogeneous catalysis art above except that the boxes or spheres used to cover the 
surface when measuring the surface's box counting dimension must be larger than 
probe gas molecules. To do this, the commercially available software was merely 
cited in the specification to teach a convenient way of producing the boxes/spheres in 
applying the box counting method, for example, to the micrographs in the 
specification (page 6, lines 8-26). 

To illustrate the box counting method, a surface is covered by spheres or boxes 
(squares) of varying sizes and the amount of spheres or boxes counted, where 
differing rules of placement of the boxes or spheres may be used as long as it is self 
consistent as illustrated in Falconer on page 40. The box count ing/fractal dimension 
is then determined merely by plotting the log of the number of boxes versus the log of 
the size of the boxes wherein the slope of the line is the box counting/fractal 
dimension in the size range selected as described by the specification on page 6, lines 
8-26 and described on page 37-40 of Falconer. The surface box counting method is 
depicted in the following illustrations 1-4. 
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Illustration 3: Idealized surface covered with larger spheres. 



Illustration 4: Log-log plot of number of spheres v. sphere size 



Upper 


sphere limit 

s 

Log sphere size 

Since applicants specifically describe and use the box counting method to 
determine the fractal/box counting dimension consistent with its ordinary meaning 
and the software is merely cited in the specification as a convenient way of producing 
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the boxes when analyzing a micrograph of a cross-section of the catalyst-substrate 
surface, Independent Claims 1, 2 and 16 and Claims dependent therefrom are definite. 

THE § 102(b) REJECTION 

Claims 1, 2 and 4-22 have been rejected under 35 U.S.C. § 103(b) as 
being anticipated by Guile et al., U.S. Pat. No. 5,716,899 (Guile herein). Applicants 
respectfully disagree and respond as follows. 

Applicants invention, in claims 1, 2 and 16, is a catalytic device that is 
comprised of a ceramic support of fused ceramic grains and a catalyst in which the 
surface is defined by a box counting dimension. The smallest box used to determine 
the box counting dimension is 1 micrometer. That is. any feature of the surface on the 
order of 1 micrometer or smaller is disregarded and as such will not impact the box 
counting dimension (i.e., will render a box counting dimension of 1 or close to 1 ). 
Claims 1, 2 and 16 make it clear that the surface is significantly fractal in nature or if 
periodic a bed of nails type structure where the surface feature is at least greater than 1 
micrometer. 

In contrast, Guile describes a bare (i.e., honeycomb without a catalyst) porous 
honeycomb that has pores having an average diameter of about 1 5 to 30 micrometers, 
(col. 4, lines 1 3-20). Guile, then, specifically describes that the pores of the bare 
honeycomb are filled with "active" material or in other words a catalyst, (col. 4, lines 
40-45). Thus, the surface of Guile now is a surface that is either dense flat 
honeycomb surface or flat catalyst surface as shown in Fig. 2a. The activated carbon 
particles have pores on the order of 1 0 to 15 angstroms, which is much smaller than 
the lower limit of the box counting dimension, (col. 5, lines 50-56). The following 
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Figures 1 and 2 illustrate the bare substrate and an idealized filled "catalyzed" 

substrate of Guile. 

Fig. 1 : Bare substrate having pores. 

Substrate surface 



Fig. 2; Guile substrate where pores are filled with catalyst particles. 


Substrate surface 



From Guile, the catalyst particles must be substantially smaller than the pores 
of the honeycomb body, otherwise they would not fill the pores, (col. 6, lines 12-14). 
The catalyst particles are ground and the median particle diameter after grinding is 
about 2-6 micrometers. It is well known that monospheres of particles closed packed 
will result in a pore diameter of less than 1 micrometer. (See page 190 of Introduction 
to the Principles of Ceramic Processing , James S. Reed, Wiley Interscience, New 
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York, NY, 1988, previously submitted). It is also well known that continuous size 
distributions of particles, such as those arising from ball milling as in Guile, will pack 
much more efficiently (the larger pores are filled by the smaller particles), id. The 
resulting Guile catalyst-substrate surface, by first principles, does not create the 
required surface of the present invention because the pores and features of the surface 
will not have features even remotely greater than 1 micrometer. That is Guile's 
catalyst-substrate surface will have pores less than 1 micrometer and as such a sphere 
or box having a dimension of 1 micrometer or larger will span such a pore and as such 
will appear as smooth surface (straight line) in the box counting method as 
exemplified in illustration 3 above. Because Guile's catalyst-substrate surface will 
appear as a straight line in the box counting method it will have a fractal/box counting 
dimension that is or is nearly 1 . Consequently, Claims 1 , 2 and 4-22 are novel. 

In addition, since Guile specifically creates a surface that fills in all of the 
pores with fine particles. Guile teaches away from the surface structure of Applicants 
invention. For this reason. Independent Claims 1, 2 and 16 and Claims dependent 
therefrom are non-obvious over Guile. 

THE § 103(a) REJECTION 

Claim 3 has been rejected under 35 U.S.C. § 103(a) as being unpatentable over 
Guile, as stated in the above paragraph (referring to the Office Action's 102(b) 
rejection), and further in view of Moyer et al., U.S. Pat. No. 5,194,154 (Moyer 
herein). Claims 23-31 have been rejected as being unpatentable over Guile in view of 
Fukui et al., U.S. Pat. No. 5,569,455. 

As stated above, Guile specifically teaches to fill in the pores with fine catalyst 
particles , which is opposite of the surface as required in independent Claim 1 , 2 and 
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16. Even if Guile was combined with Moyer, the mullite honeycomb would be filled 
with catalyst particles such that the same structure of Guile would result. That is 
Guile would fill all the pores between the mullite grains making a surface not having 
features (pores) greater than 1 micrometer and as such would not create the present 
invention. Consequently, Guile in view of Moyer fails to teach the invention, but 
teaches away from the invention of Claims 1 , 2 and 1 6. For this reason Claims 1 , 2, 
and 16 and Claims dependent therefrom are non-obvious over Guile in view of 
Moyer. 


Claims dependent therefrom are patentable. Applicants, therefore, respectfully 
request withdrawal of all rejections and allowance of Claims 1-37. 


Considering the foregoing reasons. Independent Claims 1, 2 and 16 and 



Kevin jOvilsen ^ 
Registration No. 41,510 
Phone: 989-638-6505 
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The characterization of heterogeneous catalysts or porous solids involves a description of a 
reacting interface. The physical and chemical properties of catalysts - surface area, active silc 
density and reactivity - have recently been quantified by fractal geometry. The use of fractal 
dimension rather than N 2 BET surface area provides for the determination of effective 
surface area for any adsorbate. The concept may be useful in predicting the reactivity of 
individual components in hydroprocessing. Similarly, the effects of particle size on surface 
area and reactivity, or the distribution of active sites on a catalyst, may be quantified by 
fractal geometry. Applications are discussed with specific examples. 


1. Introduction 

The geometry of a reacting interface is extremely important in most industrial 
reaction systems. In heterogeneous catalysis, the performance of a catalytic 
system is dictated by the geometry of the interface as well as by the chemical 
reactions which occur at that interface. The application of fractal geometry for 
characterizing boundaries or porous materials is currently receiving considerable 
attention. This paper discusses potential application of fractal geometry in 
heterogeneous catalysis, particularly hydroprocessing. 

The concept of spatial dimension has normally been represented by integers. In 
classical Euclidean geometry, a line, rectangle or cube are described as one, two, 
or three-dimensional. The physical world, however, involves nontrivial geome- 
tries. Mandelbrot [1] introduced the term fractal as a dimensionless number 
which quantifies the degree of irregularity or rugged n ess of a boundary. A fractal 
therefore is an object that has a fractional dimension. A rough surface, for 
example, may have an intermediate dimension between a 2-dimensional smooth 
surface and a 3-dimensional volume. The fractal dimension of a surface, D, is 
therefore a quantitative measure of the space-filling ability of that surface. 

In the field of heterogeneous catalysis, the coupling of geometry and chemistry 
introduces the phenomenon of structure sensitivity. The complex geometry of the 
interface may be described quantitatively by a fractal dimension if the interface is 
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symmetric to transformations of scale, i.e., self-similar. This simply means thai 
for a fractal, the same type of features are observed at different magnifications. 
In heterogeneous catalysis this change of scale can involve reactants of different 
molecular sues (as in the sorption of a series of adsorbates), or catalyst particlei 
of different sizes with the same reactants, or supported catalysts with varyint 
amounts of active sites. Avnir and colleagues [2,3) have introduced two simple 
equations to quantify the scaling phenomena of surfaces. If a catalyst surface is a 
fractal, the accessible surface area, A, is related to catalyst particle radius, R, by: 
A k x R (1) 

where D is the surface fractal dimension. In addition, accessible surface area is 
related to the effective cross-sectional area of the adsorbate, a m , by: 

A * lfl ~ (2) 
The surface fractal dimension is a very useful quantitative parameter for char- 
acterizing typical hydro-processing catalysts. 

Nitrogen BET specific surface area is most often used as a quantitative 
measure of surface area. The porous materials used as supports for hydrotreatin* 
catalysts - amorphous alumina and silica-alumina - have a surface fractal 
dimension between 2 and 3. The larger the fractal dimension, approaching 3 the 
more rugged the catalyst surface. In the limit, the surface will fill a volume of 
■pace. The concept of the space-filling ability of a catalyst leads to a slightly 
different approach in describing catalytic reactions. In effect, a catalytic reaction 
may not so much occur on a surface as in a surface. The active sites are normally 
viewed as elements of disorder - kinks, vacancies - when perhaps an alternate 
picture for hydroprocessing catalysts is an element of order in a disordered 
surface. 


2. Discussion 


The equations above are useful not only in measuring D but also in calculatina 
accessible surface area, effective cross-sectional area of adsorbates or even 
particle size. For reactions of large molecules, the use of accessible surface area 
for each reactant/ product is preferable to the use of N 2 BET surface area Table 
1 was calculated from eq. (2) for three hypothetical hydrotreating catalysts of 
surface fractal dimension Z> = 2.10, 2.50 or 2.90. The proportionality constant 
was determined by assuming a value of 200 m'g- 1 for the N, BET area, and then 
the accessible areas for benzene, naphthalene, and anthracene were calculated 

Similarly, by assuming an anthracene surface area of 100 m 2 g _1 the N BET 
specific surface area can be calculated as in table 2. 2 

The area occupied by an adsorbate molecule in a monolayer, a . was calcu 
lated from liquid density values [4]. The tables indicate that N 2 BET surface area 
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Table 1 

Accessible surface area 

absorbate a m (run 3 ) Calculated accessible surface area 

(mV 1 ) 

D - 2.10 D - 2.30 D = 2 90 

nitrogen 0.162 20o" 200 200 " 

benzene 0.305 194 171 150 

naphthalene 0.383 192 176 136 

anthracene 0.410 191 159 132 


is a poor parameter for characterizing the three catalyst surfaces. Considering a 
complex mixture of reactants (such as hydrogen, benzene, napthalene and 
anthracene) it is possible to choose a catalyst that has relatively more surface 
accessible to anthracene, for example, by choosing the catalyst with the lowest 
surface fractal dimension. The underlying assumption, of course, is that the active 
sites are evenly distributed throughout the catalyst surface. 

For catalytic hydrocracking then, the catalyst with the lowest surface fractal 
dimension (assuming equivalent N 2 BET areas) should present the highest 
accessible area for asphaltenes. Hydroprocessing catalysts are often described as 
sponges or as porous solids with very high "internal" surface as opposed to 
"external" surface. High "internal" surface is ascribed to a pore structure where 
the classical mathematical model of a pore is a cylinder with smooth walls This 
model may be described by a Menger sponge model [1,5]. The alternate picture of 
these materials is as a very rugged surface which is space-filling to a degree 
described quantitatively by a surface fractal dimension. The relationship between 
surface fractal dimension, Z>, and average cylinder-equivalent pore diameter has 
been explored [6]. The surface fractal dimension, Z>, is clearly a useful quantita- 
tive parameter for surface characterization. 

An interesting catalyst has been described for the ASVAHL hydrodemetalliza- 
tion (HDM) process [7]. This is a burr-shaped catalyst with no "internaT surface 
as opposed to a sponge. The branching of this burr appears to be self-similar over 
a range of scale. The outside of the burr would appear smooth to a large 
asphaltene molecule. The initial reactions of HDM and cracking occur at the 
outside of the catalyst particle and it is here that metals are deposited. Smaller 


Table 2 


adsorb ate 

Calculated N 3 

BET surface area 







Z>*2.10 

D - 2.50 

D ~ 2.90 

anthracene 

100 

100 

100 

nitrogen 

105 

126 

152 
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molecules penetrate the burr where further reactions may occur. The entire 
catalyst then would appear to open up to expose new surface to geometrically 
smaller molecules. In describing this catalyst as a porous material, one would 
model it as a solid with pore sires ranging from very large at the outside edges to 
small in the center of the bum 

Equation (1) quantifies the effect of particle size on apparent surface area. For 
a high D particle, the apparent surface area is independent of R. For particles of 
intermediate surface fractal dimension (such as most transition metal oxides), this 
equation gains importance as it does in fine particle science (8J. Nitrogen BET 
data for microporous materials should be treated with caution. Applications of 
fractal geometry in adsorption studies on microporous petroleum coke have been 
examined [8b]. 


3. Reaction dimension 

The numerous parameters which dictate the chemical reactivity of a fractal 
boundary may be described by an empirical scaling exponent, D R the reactive 
surface fractal dimension. This can be interpreted as the effective fractal dimen- 
sion of an interface towards reaction. It has been demonstrated [9] that the initial 
rate of reaction, u f with the surface of a particle i.e., thermal decomposition or 
gasification is related to particle size by: 

= (3 ) 

The reactive surface fractal dimension is a measure of the activity of the reactive 
surface which is a sub-set of the total particle surface. Farin and Avnir have 
reviewed several examples of eq. (3) and discussed the relationship between D R 
and D [9]. Equation (3) should also find application in heterogeneous catalytic 
reactions. 

The catalytic activity of metal catalysts dispersed on high surface area supports 
has been shown [10] to exhibit similar scaling behaviour; 

a = k 4 R$* (mol s" 1 per particle) (4) 
a = k i Rp M ~ 2 (mo\ s" 1 per exposed surface site) (5) 
a^M^-^rnols- 1 g" 1 ). (6) 

In eqs. (4) to (6), R p represents the radius of the metal particles on the support 
and D K is the reaction dimension. The reaction dimension provides a quantita- 
tive measure of structure sensitivity. Farin and Avnir [101 have noted a range of 
values from D R = 0.4 to D n - 6. The reaction dimension is affected not only by 
specific molecule/metal interactions but also by the support morphology. 

When catalytic activity per exposed surface atom is independent of metal 
crystallite size, the reaction is structure insensitive. This occurs when D R = 2, eq. 
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(5). In terms of activity per particle, cq. (4), the activity scales with R 2 . Therefore, 
all surface atoms contribute equally to catalytic activity regardless of crystal 
plane. The reaction dimension is easily determined from log-log plots of eqs. (4), 
(5) or (6). The real value of Z> R , however, occurs when specific active site 
populations can be assigned to Z> R values. The easiest specific sites to visualize 
may be corners, edges, and planes. Where D K is less than 2 and the surface of the 
metal crystallite has dimension 2, a subset of active sites, exists which grows 
slower than the total number of surface sites, n s . The opposite case. Z> R greater 
than 2, means that the subset of active sites grows faster than the total number of 
surface sites. The performance of this particle size scaling analysis indicates that 
if the particle is magnified by a factor of N t then a total of N° features appear 
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4. Applications of reaction dimension 

A great deal of published literature can be examined to demonstrate a scaling 
analysis of typical Mo/A! 2 0 3 catalysts. The rate of hydrodesulphurization (HDS), 
for example, may be related to the effects of molybdenum loading. Where eq. (6) 
relates catalytic activity to metal crystallite size in supported catalysts, it is 
perhaps easier to relate HDS activity to molybdenum loading for hydrotreating 
catalysts. The data of Okamoto ct al. [11] are plotted in fig. 1 for thiophene HDS 
This is a log-log plot of thiophene conversion (400 °C, atmospheric pressure) 
against molybdenum loading. The loading has been expressed as the number of 
molybdenum atoms per square nanometer of N 2 BET surface area. This unit of 
loading was chosen as a better approximation of surface site density to facilitate 
comparisons of different alumina supported catalysts. Figure 1 exhibits regions of 
structure sensititivity for thiophene HDS above and below a molybdenum load- 
ing of 2.7 Mo nm 2 . The reaction dimensions for the two linear regions 
calculated from the slope of the log-log plot using eq. (6) are 5.2 and 3.7 for the 
low and high loadings. 

Similarly, the data of Thomas et al. [12) for thiophene HDS (400 °C, atmo- 
spheric) are presented in fig. 2. These data were generated to observe the effects 
of impregnation methods (i.e., wet and dry) on catalyst activity. Clearly, when 
plotted in terms of m* thiopbHhe converted per gram-catalyst second, the various 
preparation techniques produce final catalysts with the same structure sensitivity. 
The reaction dimension, £> Rf determined from these data is 4.5 below a loading 
of 2.7 Mo nm 2 and 3.4 above: Bachelier et al. [13] have examined the effect of 
molybdenum loading on thiophene HDS (420 °C, atmospheric) and on oxygen 
sorption at 60 °C. Catalysts were prepared by a normal aqueous impregnation of 
gamma alumina using ammonium heptamolybdate. The data are shown in fig. 3 
as a log-log plot of reaction rate (mol g^h" 1 ) or oxygen adsorption (umol g" 1 ) 
against loading. Two linear regions were again observed. For thiophene HDS, a 
reaction dimension of 5.2 is calculated for loadings below 2 Mo nm" 2 and 3.7at 
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Mo leodtng (Me / nm*) 

Fig. 1. The conversion of ihiophene over a series of Mo/Al ,0, catalysts; Q, ref. [11a]; + , ret (lib]. 


higher loadings. The oxygen adsorption data followed similar trends with calcu- 
lated reaction dimensions of 5.0 and 3.2 for the low and high loadings respec- 
tively. ' y 

The three examples demonstrate that the structure sensitivity of thiophene 
HDS can be quantified by a reaction dimension, £> R . In all cases, there is a 
change in structure sensitivity at molybdenum loading between 2 and 3 Mo 
nm It is anticipated that this point coincides with the beginning of molybdenum 
aggregation with increased loading. 

It is apparent [14-16J that MoO, can be detected by laser Raman Spectroscopy 
(LRS) only at loadings above 3 Mo nm" 2 after calcination at 500 °C. Extended 
X-ray adsorption fine structure (EXAFS) indicated no evidence of MoO, in a 
catalyst of 2.4 Mo nm" 2 loading [17.18]. The EXAFS study suggest a large degree 
of disorder with very small amounts of crystalline material, i.e., less than 1 nm 
diameter. At loadings below which Mo0 3 or Alj(Mo0 4 ) 3 are detectable by LRS 
or EXAFS, an interaction species has been described [14,19]. A tctrahedral 
coordination for molybdenum predominates at low loadings, tetrahedral and 
octahedral coordinations at 2.8 Mo nm -1 . An increase in both octahedral and 
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Mo loodlrxj (Mo / nm°) 
Fig. 2. The rale or thiophene HDS for a series of Mo/Ai 2 0 3 catalysts; Qdry impregnation- 


impregnation: H. effect of pH. B; boehmite; ref. (12J. 


5.0 


, wet 


tetrahedral molybdenum with increased loading was observed [14J. A maximum 
in tetrahedra] coordination occurred at 1.7 Mo nm" J and in octahedral at 4 5 Mo 
nm loadings. Tetrahedra! molybdenum was thought to occur in isolation i e 
well-dispersed whereas octahedral molybdenum was thought to exist in associa' 
tion with other Mo groups. 

The data presented in the figures clearly describe an empirical scaling phenom- 
enon. The application of a fractal descriptor, £> R , results in a quantitative 
description of catalyst structure sensitivity. This is useful in comparing catalysts 
prepared by different methods using various alumina supports. In addition there 
is evidence to suggest that the molybdenum loading may be related to molybdenum 
dispersion, aggregation, or the size of ordered regions on the catalyst surface 
There may be examples in the literature where rates of reaction for a series of 
reactants are presented. In some cases it has been stated that large molecules 
react at a lower rate than small molecules due lo steric hinderance. The reaction 
rates are most often presented as moles per gram or moles per m 2 N, BET 
surface area. As indicated in lables 1 and 2 however, large molecules simply do 
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Fig. 3. Thiophcnc HDS activity, a and ox>gcn chcmisorption capacity. + for a srries »r 

not "sec" the same quantity or surface as small molecules. Where rate is 
controlled by molecular/ surface geometry, the initial rale of reaction is con 
trolled by the effective reactant cross-sectional area, a : 


-D H /2 


(7) 


In this example, a log-log plot of rate of reaction against reactant cross-sectional 
area would give a straight line of slope - 

Wilson et al. [20J present data for the hydrogenation of a synthetic crude 
middle distillate over a commercial Mo-Ni/Al 2 0 3 hydrotreating catalyst The 
feed and products were accurately characterized in terms of alkylbenzene benzo 
cycloparaffin and benzodicycloparaffin fractions. The rate of hydrogenation 
obviously decreased with increased reactant size. When plotted according to eq 
(7), the data indicate a reasonable relationship, as shown in fig. 4. The effective 
cross-sectional areas of fractions were approximated as before [4] using benzene 
iciralin. and octanthrene as model reactants. In spite of the approximations in the 
rate measurements and estimates for a m9 the trend is apparent. In this discussion 
the scaling parameter is the size of the reactant molecule. The measured reaction 
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1 " r 

3.047E-19 4.051 E-1 9 4.905E-19 

Effectlvo area ( m l ) 

Fig. 4. The rate of hydrogenalion as a function or reactant size for hydroproccssing a svmheiic 
crude distillate over Mt>-Ni/AI ,0,: a 380°C: + . 360°C; 0. 340°C; ref. [201. 

dimension, D K = 2, suggests that the reaction is not structure-sensitive with 
respect to reactant size. The catalyst, therefore, is equally accessible to all three 
molecules for hydrogenalion. Since the accessible surface area varies with ad 
sorbate size for a porous catalyst, eq. (2). the conclusion is that catalyst reactivity 
towards hydrogenetion occurs only on the external parts of the catalyst Similar 
scaling analyses for other reactions such as hydrodesulphurization, hydrodc- 
nitrogenation, hydrodemetallization or asphaltene conversion would provide use- 
ful information on the structural sensitivity of these important reactions. 


5. Conclusions 


The properties of typical hydroproccssing catalysts have been discussed from 
the standpoint of fractal geometry The potential of the fractal approach for 
characterizing catalysts and catalytic reactions has been demonstrated The 
principal conclusion is that hydroprocessing catalysts are more accurately quanti- 
fied by surface fractal dimension, £>. rather than by N 2 BET. Equations relating 
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the apparent surface area to particle size and adsorbate cross-sectional area are 
also appropriate for characterizing hydrotreating catalysts. The concept of a 
reaction dimension appears to be a promising parameter for characterizing the 
reactivity of heterogeneous catalytic systems. 
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Abstract 

To investigate the effect of surface roughness on diffusion and reaction in porous 'amorphous' catalysts, a realistic 
model of roughness is necessary. By combining small-angle and wide-angle X-ray scattering, the suggestion from 
adsorption experiments that the surface is often fractal can be confirmed, since these techniques allow to measure the 
surface and bulk properties over a wider range of scales, so that the complete scaling regime, if present, can be 
revealed. The fractal dimension can be determined, as well as the inner and outer cutoff of the fractal scaling regime. 
No significant fractal surface roughness is found for zeolites, but is shown for such materials as a mesoporous 
Y-Al 2 Oj, catalytic reforming and hydrotreating catalysts. Results from the simulation of aggregation phenomena 
similar to those occurring during the sol-gel synthesis used in catalyst preparation lead to additional insists. The 
fractal surface roughness of a catalyst or catalyst support on molecular scales has a significant influence on the 
conversions and selectivities of chemical reactions, especially when they are limited by Knudsen diffusion. Tins effect 
is felt up to the scale of industrial reactors, which implies that changing the fractal surface roughness can improve the 
performance of an industrial process. © 2001 Elsevier Science B.V. All rights reserved. 

Keywords: Small-angle X-ray scattering; Fractal; Porous catalyst; Diffusion; Reaction 


1. Introduction 

Recent years have seen an increased use of 
zeolites as catalysts, catalyst supports, or for sepa- 
rations [1-4]. Mesoporous sieves such as MCM- 
41 are extending the range of controlled pore 
structures into the mesopore regime and beyond 
[5,6]. Nano-controlling reactions through an im- 
posed pore network architecture is one of the 
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main reasons driving this research. Nevertheless, 
cheap disordered amorphous mesoporous cata- 
lysts and catalyst supports are still most fre- 
quently used in the chemical industry [2-4J. They 
are also not as 'amorphous' as their etymology 
suggests, since their surface or pore volume may 
contain an internal symmetry as well. Characteriz- 
ing this internal symmetry — self-similarity 

and studying its influence in catalysis, is the object 
of this paper. 

Zeolites and M41S periodic mesoporous silicas 
obey a translational symmetry, i.e. there are re- 
peating units in one, two or three dimensions. 
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Pfeifer and Avnir's pioneering work in the eighties 
[7-9J demonstrated the possibility for the surface 
of many 'amorphous' porous materials to obey 
another internal symmetry, that had been pro- 
posed earlier by Mandelbrot [10] as a remarkably 
universal symmetry in nature, that of fractal self- 
similarity. Through the adsorption of molecules 
of different sizes, acting as probes of surface 
roughness, Avnir and Pfeifer [7-9] showed that in 
many cases, the internal surface area S of porous 
materials is size-dependent in a particular way: it 
is a power law of the molecular probe size S y i.e. 
S~6 2 ~ D , where D lies between 2 and 3, both 
limits included. This was interpreted as being a 
result of the tortuous nature of the surface, very 
similar to a natural coastline, as Adamson [11] 
already intuitively noted: large fjords or irregular- 
ities along the surface are bordered by smaller 
ones, and even smaller fjords border the latter in 
a similar way, and so on, from a largest scale (the 
outer cutoff) to a smallest (the inner cutoff) [7-9]. 
This is schematically shown in Fig. 1 . This invari- 
ance under magnification is typical for fractals 
[10]. The exponent Z), the fractal dimension, is a 
measure of the space-filling capacity of the sur- 
face; it is 2 for a Euclidean surface, such as the 
smooth surface of graphite, and it can be as high 
as 3 for the space-filling surface of some silica or 
alumina gels. Mandelbrot's [10] self-similarity of 
coastlines and other meso- and macroscopic 
shapes was thus extended into the molecular 
range. 


The results of Pfeifer and Avnir inply that, if 
the surface is fractal on molecular scales, the 
surface area measured by nitrogen ad»rption, the 
most commonly used technique, is insufficient, in 
view of the fact that in many chemical reactions 
or sorption processes using these sane porous 
materials, the species have a very different size 
To which extent it is insufficient can hr seen from 
simulation results, as presented is the next 
section. 

A lot of work has since concentrated on the 
characterization of molecular fractal surface 
roughness of porous materials [12-16]. One of the 
most frequently applied techniques has been ad- 
sorption measurements using differer* probes. It 
was quickly recognized, however, that this tech- 
nique has its limitations, as is shortly reviewed in 
the next section. However, alternatives bypassing 
these limitations exist, one of them being Nei- 
mark's highly practical thermodynamic method 
[17-19], which can be applied to a conventional 
nitrogen adsorption/desorption experiment. An- 
other powerful alternative is to combine small-an- 
gle and wide-angle X-ray scattering 
(SAXS/WAXS), which is the technique used in 
this paper. The internal surface roughness on 
molecular scales can also be probed byNMR [20], 
and using multilayer adsorption [21-231, amongst 
others [12]. Because of the less straightforward 
interpretation of the results, the latter techniques 
are preferably used in conjunction with one of the 
former. 


Fig. I. Schematic representation of monolayer adsorption on the fractal surface of a catalyst pore. The small molecules \ 
have access to a larger surface than the large molecules on the right. 
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A recent paper by the author reviews the effect 
of surface roughness on diffusion and reaction in 
porous catalysts [15]. Some of the conclusions will 
be repeated here, since they show the potential of 
fractal models to not only perform reaction engi- 
neering calculations, but also, by tuning the geo- 
metrical properties of the surface, to control 
conversions and selectivities of chemical pro- 
cesses. First, attention will be paid to the fractal 
surface characterization of porous catalysts using 
SAXS/WAXS. 


2. Characterization of the surface roughness of 
porous catalysts 

To investigate the geometrical heterogeneity of 
the surface of porous catalysts on molecular 
scales, several techniques are available, the princi- 
pal of which were mentioned in the introduction. 
Monolayer adsorption of molecules with a differ- 
ent size S is an elegant technique that has first 
indicated the possible fractality of the rough sur- 
face of many porous catalysts. The number of 
molecules in a monolayer on a fractally rough 
self-similar surface scales like N ~ 6~ D , so that 
the surface area should scale like S~ d 2 ' D for a 
self-similar surface with fractal (adsorption) di- 
mension D. However, there are some disadvan- 
tages to this method as well [15-17J. First, it is 
not easy to carry out a series of experiments with 
molecular probes with a sufficiently differing size, 
so that power law scaling can only be checked 
within a limited range. The use of very large 
molecules, such as polymers, could extend the 
range of investigation, but is often cumbersome 
and the question is shifted to a correct representa- 
tion of the conformation of the adsorbed 
molecules. Checking power law scaling within a 
narrow range of data, say less than a decade, is 
dangerous; values different from 2 or 3 may be 
too quickly interpreted as being a result of surface 
irregularity, since the statistics are not sufficiently 
accurate. In addition, if the surface is indeed 
fractal, its scaling regime can rarely be obtained. 
Finding the scaling regime is useful for practical 
applications in, e.g. catalysis, and to help inter- 
preting the results. A second problem with ad- 


sorption measurements is the possibility of 
different chemical interactions between different 
probe molecules and the surface. In that case, fee 
results are a reflection of chemical, rather than 
geometrical heterogeneity, or a mixture of both 
Separating the effects of different types of hetero- 
geneities is difficult [24,25]. 

Nevertheless, adsorption methods using a singfe 
probe, typically nitrogen, can be used to investi- 
gate surface roughness. Pfeifer and Avnir first 
proposed to compare adsorption measurements 
carried out on particles of different sizes, to ex- 
plore how the specific surface area, measured 
with, say, nitrogen, scales with particle size R. A 
fractal self-similar surface with dimension D 
would lead to the scaling law S~R°-*. How- 
ever, this result is only valid if the particles aic 
formed through some self-similar process leading 
to an outer cutoff that is proportional to R as 
well. Neimark introduced an easy to use method 
based on capillary condensation of the inert nitro- 
gen on a fractal surface, for which S^al'*, 
where a c is the mean radius of curvature. The 
latter is related to the relative pressure P/P Q 
through the Kelvin equation, while S at any rela- 
tive pressure can be found by integration of the 
isotherm from the current value to the maximum 
value for P-+P 0 according to the Kiselev equa- 
tion [17-19]. In fact, the method is quite similar 
to porosiinetry, since it is applied to adsorption m 
the pore-filling regime. This thermodynamic 
method can be used for mesopore irregularities 
larger than 1-2 nm. 

An alternative method to adsorption, as al* 
ready indicated in the introduction, is small-angle 
X-ray scattering (SAXS). This is a useful alterna- 
tive that can provide interesting information over 
a wide range of length scales, and which is a 
purely physical technique, not prone to network 
hysteresis effects or energetic effects that may 
affect adsorption methods. For an introduction to 
SAXS, see [26-28]. The intensity / of an X-ray 
beam, scattered by the sample at the same wave- 
length X as the source, is measured as a function 
of the scattering angle 6. The scattering is a result 
of the position dependence of the electron density. 
The intensity / measured as a function of the 
length of the scattered wave vector, q = 4n/A 
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sin(0/2), yields information about structural (or, 
actually, electron density) correlations on scales 
5 ~ l/q (viz. Bragg's law). Because the angles at 
which the measurements are performed with a 
SAXS apparatus are much smaller than with 
XRD (down to a few minutes), information on 
colloidal scales or scales between about 1 and 100 
nm (sometimes even larger) can typically be ob- 
tained. Crossovers and features of a specific size 
in this range can be measured with SAXS, which 
is, therefore, popular in polymer science, and for 
the measurement of the size of micelles and col- 
loidal particles. It has also been used for the 
determination of particle size (distributions) in 
catalysis, by filling the pores with an electron 
density matching liquid, to avoid scattering at the 
pore surface [29,30]. However, also pore size dis- 
tributions and surface roughness can be conve- 
niently probed. In this context, it is useful to 
remember that fractality is a symmetry, leading to 
specific power-law correlations. Through scatter- 
ing techniques, these correlations are measured in 
the Fourier domain, where they also lead to 
power laws. A nice overview of the application of 
small-angle scattering to fractal morphological in- 
vestigations is given by Schmidt [31,32]. For a 
self-similar surface fractal, Bale and Schmidt [33] 
showed that f~q D ~ 6 , which is an extension of 
Porod's I~*q~ 4 law for smooth surfaces. To 
derive this equation, they used a two-phase ap- 
proximation: the porous material consists of two 
phases with a constant electron density, and the 
interface is a self-similar fractal. In addition, the 
solid is assumed to be representable by a set of 
rigid solid spheres. Apart from a vast literature on 
fractal aggregates [34,35], several solid materials 
were tested using SAXS (and SANS, a similar 
technique using neutrons), such as porous rocks 
(14), Vycor porous and controlled-pore glass [36], 
as well as silica xerogels and coal particles [37], to 
give a few examples. The technique seems equally 
applicable to the industrial mesoporous catalysts, 
the problem of interest to this paper. 

One potential drawback of scattering methods 
is that the investigated surface of a porous mate- 
rial is not necessarily the surface accessible from 
the outside, i.e. the surface accessible by diffusing, 
adsorbing and/or reacting molecules. Part of the 


pore surface may indeed be occluded. The way in 
which the pores and their surface are formed 
should indicate whether inaccessible pores consti- 
tute an important part of the material and 
whether the inaccessible and accessible surface 
have different morphologies. In case of doubt, 
density matching liquid pore filling can be used to 
separately study the occluded pores, a technique 
especially useful to study porous carbons. This 
was not done for the catalysts studied here, be- 
cause their sol-gel based synthesis should not 
lead to a significant occluded porosity or to a 
different surface morphology of occluded and ex- 
ternally accessible surface; the reaction studies in 
the next section can, therefore, make use of the 
(fractal) surface properties derived in this section. 

The focus here is on a number of mesoporous 
alumina supported catalysts, used in refinery and 
petrochemical processes. Results for a PtRe/AIjOj 
catalyst used in the catalytic reforming of naphtha 
are shown in Fig. 2 by way of example. This 
figure combines results obtained with three differ- 
ent apparatus, to cover a range of scales going 
from Angstroms to about 50 nm. Results for the 
largest scales were obtained at very small angles, 
using synchrotron radiation at the Institut Laue- 
Langevin in Grenoble — the downward deviation 
at the smallest angles is a result of the beam stop, 
which is blocking the (overpowering) radiation 
passing straight through the material. The data at 
average scales were obtained using X-rays from 
the synchrotron in Daresbury. Finally, the data at 
the smallest scales, corresponding to the largest 
angles, cannot be obtained by SAXS, but were 
acquired using wide-angle X-ray scattering 
(WAXS) in the Shell Research and Technology 
Center in Amsterdam. A nice mutual overlap of 
the results can be observed, covering an impres- 
sive 6-7 orders of magnitude in intensity. By 
using WAXS, the for catalysis so interesting 
molecular range can even be probed. A clear 
power-law scaling is seen in a relatively broad 
regime; the inner and outer cutoff of the scaling 
regime can easily be calculated from the knees in 
the curve, d^^l/q^ and S m3x ^n/q min [3I|. 
The smallest scale is indeed seen to be of the order 
of very small molecules, a few Angstroms in size, 
while the largest scale is similar to the average 
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Fig. 2. SAXS/WAXS measurement results for a PtRe/Al 2 0 3 catalytic reforming catalyst. 


mesopore size of this material and the size of the 
aggregates it consists of, ~ 6 nm. The results also 
confirm those obtained for another reforming cat- 
alyst within a narrower regime, using a Kratky 
camera with a Cu-K^ X-ray source. A pretreated, 
presulphided catalyst that has seen the reactor 
feed, gives similar results. The Pt and Re content 
of the catalyst are 0.3 wt.% each, and it can be 
verified that such a small content has very little 
effect on the results (note that the particle size can 
be found using electron density matching, as men- 
tioned earlier [29,30]). All this points toward a 
consistent measurement of the geometrical surface 
characteristics of the catalyst. 

Surprising was the absolute value of the slope 
in the scaling regime, however, which was found 
to be 4.33 with high accuracy. This is significantly 
larger than 4, corresponding to a smooth surface 
(Porod's law), and represents a very steep de- 
crease of the intensity with q~l/d. A value 
smaller than 3 would have corresponded to a 
mass or volume fractal (I^q~°) [31,32], while a 
value between 3 and 4 is traditionally attributed 
to a fractal surface, using Bale and Schmidt's 
expression I~»q D - 6 . Values larger than 4 have 
rarely been reported. Experiments with CoMo/ 
Al 2 0, and NiMoP/AI 2 0 3 hydrotreating catalysts 
led to slightly different values of the slope and 
cutoffs, but a very similar behavior with again a 


slope larger than 4. An alumina support led to an 
even higher slope of 4.6, again using mutually 
corroborating, reproducible SAXS and WAXS 
results. 

Fig. 3 shows the scattering results for a Pt-ul- 
trastable Y faujasite used for hydrocracking of 
gasoil. Zeolites like faujasite are crystalline alu- 
mino-silicates, in contrast to the amorphous 
mesopores materials investigated earlier. The dea- 
lumination of faujasite by steaming leads to the 
formation of mesopores. Their size of around 3 
nm (knee at q - 1 nm '~ 1 in the SAXS/WAXS) 
nicely coincides with a typical value mentioned in 
the literature for cages in this catalyst [38| The 
slope for q>\ nm" 1 is 4.0, consistent with 
Porod's law for non-fractal surfaces. 

To find an explanation for the origin of the 
power larger than 4 for the mesoporous aluminas, 
the origin of Bale and Schmidt's scaling law 
should be recalled. For such materials, the two- 
phase model could be an oversimplification. A 
representation of the alumina as consisting of a 
collection of spherical aggregates with fractal sur- 
face roughness up to the order of the size of the 
aggregates themselves (as observed) implies a 
gradual power-law decay, rather than a sudden 
drop of the density from the center of the aggre- 
gate to its outside. In this context, it is interesting 
to mention that scaling with a power larger than 
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4 was also found by Schaefer [39J for surfaces 
covered with a tensio-active compound, and by 
Schmidt et al. [40] for reversed phase silicas with 
grafted alkane layers that have a density dropping 
as a power law with the distance to the silica 
surface. A more refined analysis again leads to a 
power larger than 4. Moreover, based on a solid- 
on-solid model of surface roughness, a detailed 
derivation by Wong [41] led to a law / ~ q 2 + ° for 
certain fractally rough porous media, correspond- 
ing to powers between 4 and 5, as also seen in the 
experiments presented here. Further investigation 
of this interesting fractal scaling is underway. 

It should be stressed that the fractal characteri- 
zation should always, whenever possible, be done 
with the source or preparation of the material in 
mind. Many mesoporous amorphous catalysts 
and supports are prepared through a series of 
steps that starts with a sol-gel synthesis. Aggrega- 
tion in solution, far from equilibrium, often leads 
to fractal structures or particles with a fractal 
surface. This is shown both from experiments and 
simulations [34]. Subsequent filtering, drying and 
calcination may modify the texture and surface 
structure, and generally densify the particles. Also 
electron microscopy images frequently reveal a 
rough, disordered structure. A fractal surface 
within a finite scaling regime, corresponding to 


the sizes of the elementary particles and the aggre- 
gates, is, therefore, a good possibility, which may 
then be tested by the described methods, 


3. Effect on diffusion and reaction 

In chemical reactions catalyzed by a porous 
(supported) catalyst, molecules diffuse through 
the pore space, collide with its walk, and react on 
the surface [2,4,42]. Molecular scale fractal sur- 
face roughness of the internal catalyst surface 
may therefore, have an influence on diffusion and 
reaction phenomena. 

First, it influences the surface area available to 
reactions occurring on the fractal surface [12]. The 
effect of a fractal surface morphology on molecu- 
lar diffusion toward this surface, and reaction on 
it, has been studied by several authors [43-47]. A 
more detailed list of references is included in 
[45,48,15]. 

Second, it influences the Knudsen isffusivity, 
frequently the dominating diffusion mechanism 
for gas phase reactions in mesoporous catalysts 
[42,49]; since many reactions of industrial rele- 
vance are diffusion limited, this influence trans- 
lates into an effect on conversions and selectivities 
[42). Knudsen diffusion is the dominating diffu- 
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Fig. 3. SAXS/WAXS measurement results for an ultrastable faujasite Y, used for hydrocracking of gasoil. 
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sion mechanism when the mean free path of the 
gas molecules is of the same order of magnitude 
or larger than the pore size. In that case, collisions 
of the molecules with the walls are more frequent 
than intermolecular collisions. The difTusivity can 
thus be expected to be influenced by the surface 
morphology. A light molecule diffuses more 
quickly than a heavy molecule, but is usually also 
smaller. If the roughness is minor compared with 
the pore size, the size effect will be small (50,51), 
yet, if it is important, the collision frequency of 
the molecules with the walls of the irregularities 
will be significantly higher for small molecules 
that have access to more and smaller fjords than 
large molecules [15,52,53]. A mean-field approxi- 
mation shows that the Knudsen difTusivity D K is 
approximately inversely proportional to the acces- 
sible surface area [52], so that D K (S) ~ D K0 S fD ~ 2 
[there is a misprint in [15], equation 7], where D K0 
is the Knudsen difTusivity in a smooth pore and d' 
is the effective diameter of the molecules 8 divided 
by the outer cutoff of the fractal scaling regime, 
as measured by, e.g. SAXS. More detailed calcu- 
lations, based on a first-passage time calculation, 
lead to a two-parameter equation D K (d)~~D KO f 
[1 + a(l — S'*)] § where a and ft can be expressed 
as a function of the fractal dimension of the 
surface, Z), and a second parameter, the 'return 
probability', which depends on the shape of the 
fjords or irregularities [53,54,15]. 

Because of the dependence of the Knudsen 
diffusivities and the accessible surface area on 
molecular size, fractal surface roughness has an 
important effect on diffusion limited reactions. 
Models to calculate this influence are presented in 
[52,55,56] and summarized in [15]. The signifi- 
cance of the fractal roughness effect can be 
demonstrated for the case of an industrial four-re- 
actor unit for the catalytic reforming of naphtha, 
a process used to increase the octane number of 
gasoline. The calculations are based on a reaction 
network that consists of 86 reactions between 29 
component lumps, and includes a detailed de- 
scription of multi-component diffusion inside the 
porous catalyst pellets, using the Stefan-Maxwell 
based dusty gas equations [49,42] with all binary 
molecular diffusivities and fractal roughness cor- 
rected Knudsen diffusivities. Simulations of diffu- 
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Fig. 4. Octane number of the gasoline leaving an industrial 
catalytic reformer, as a function of the fractal (intension of the 
internal surface of the PtRe/AI 2 0 3 catalyst. The hydrogen 
yield at the outlet of the reformer is shown on tfte same graph. 
The results indicate that a fun her increase of the octane 
number at very high fractal dimensions is accompanied by 
excessive, undesirable hydrocracking and hydragenolysb. 

sion and reaction inside the pellets arc performed 
all along the axis of the four fixed bed reactors in 
series. The octane number of the gasoline in- 
creases thanks to aromatizations and branching 
isomenzations. Only the ultimate result of the 
detailed calculations is given in Fig. for more 
details, the reader is referred to [55,56). Because of 
the higher accessible surface at higher fractal di- 
mensions of the alumina surface, which contains 
acid sites catalyzing some of the reactions, the 
octane number steadily increases with increasing 
aromatization of the feed. However, the hydrogen 
yield peaks around D = 2.6, as a result of the too 
significant effect of undesired side-reactions in- 
volving hydrogen (hydrocracking and hy- 
drogenolysis) when the fractal dimension of the 
internal catalyst surface is too high. The reason 
for this is the relatively higher accessibility of the 
rough surface to the small hydrogen molecules, so 
that reactions converting hydrogen become im- 
portant at high fractal dimensions. The effect 
shows that a change in the catalyst surface mor- 
phology through a proper modification of the 
porous alumina synthesis may lead to an optimal 
increase in the octane number of the gasoline. In 
other words, optimization is possible by tuning 
the surface morphology. 
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4. Conclusions 

It was shown how a combination of small-angle 
and wide-angle X-ray scattering, as an extension 
and alternative to adsorption studies, can be used 
to investigate the fractal surface morphology of 
industrial porous, amorphous catalysts used, e.g. 
in catalytic reforming and hydrotreating. This 
fractal roughness is in line with, but cannot di- 
rectly be inferred from, the preparation of these 
materials. Power law scaling with a slope larger 
than 4 was discovered and a tentative explanation 
was given. The results allow to derive the com- 
plete fractal scaling regime, including inner and 
outer cutoffs, which is not usually possible from 
adsorption experiments alone. Another advantage 
of scattering techniques is their physical basis, 
independent from specific interactions between 
molecules and the surface or pore volume. The 
effect of fractal surface roughness on diffusion 
and reaction phenomena can be important, as was 
briefly shown for the example of catalytic reform- 
ing of naphtha, for which an optimal increase of 
the octane number and of the hydrogen yield can 
be predicted based on a change of the alumina 
fractal surface morphology. 
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